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Abstract 

Available throughput performance studies for mobile ad hoc networks (MA¬ 
NETs) suffer from two major limitations: they mainly focus on the scaling 
law study of throughput, while the exact throughput of such networks re¬ 
mains largely unknown; they usually consider the inhnite buffer scenarios, 
which are not applicable to the practical networks with limited buffer. As 
a step to address these limitations, this paper develops a general frame¬ 
work for the exact throughput capacity study of a class of buffer-limited 
MANETs with the two-hop relay. We hrst provide analysis to reveal how 
the throughput capacity of such a MANET is determined by its relay-buffer 
blocking probability (RBP). Based on the Embedded Markov Chain Theory 
and Queuing Theory, a novel theoretical framework is then developed to en¬ 
able the RBP and closed-form expression for exact throughput capacity to 
be derived. We further conduct case studies under two typical transmission 
scheduling schemes to illustrate the applicability of our framework and to 
explore the corresponding capacity optimization as well as capacity scaling 
law. Finally, extensive simulation and numerical results are provided to vali¬ 
date the efficiency of our framework and to show the impacts brought by the 
buffer constraint. 
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1. Introduction 

The mobile ad hoc network (MANET) represents a kind of self-organizing 
network architecture, which consists of mobile nodes communicating with 
each other without centralized infrastructure and management [1]. Since 
MANET can be deployed and reconhgured rapidly at very low cost, it serves 
as an appealing candidate for many critical applications, such as disaster 
relief, battleheld communication and emergency rescue [2], To support the 
design and applications of MANETs, the studies on their fundamental per¬ 
formance have been extensively reported. Despite much research activity, 
however, the lack of a general capacity theory for MANETs is still a long 
standing open problem and becomes an obstacle on the development and 
commercialization of such networks [3, 4], 

Since the pioneer work of Gupta and Kumar [5], extensive scaling law 
(i.e. order sense) results on the capacity of ad hoc networks have been re¬ 
ported in literature, which mainly focus on the study of asymptotic per node 
throughput behavior as the number of network nodes increases. The results 
in [5] indicate that for a static ad hoc network, its per node throughput di¬ 
minishes to zero as the number of network nodes tends to inhnity. Later, 
Grossglauser and Tse demonstrated in [6] that with the help of node mo¬ 
bility, a 0(1)^ constant per node throughput is achievable in a MANET. 
Inspired by the seminar work of [6], lots of studies have been devoted to the 
analysis on the scaling laws of MANETs throughput under various mobility 
models [8, 9, 10, 11, 12], and under various network scenarios [13, 14, 15, 16]. 
For a survey on the scaling law results of MANETs throughput, readers are 
referred to [17] and references therein. Although scaling laws are helpful to 
grasp the general trend of MANET performance, they provide a little insight 
into its real achievable throughput. In practice, however, a thorough under¬ 
standing on the real achievable throughput of MANETs is of great concern 
for network engineers, since it serves as an instruction guideline for network 
design and optimization. 
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Please kindly refer to [7] for the notations of order results. 
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By now, some initial and helpful results are available on the exact through¬ 
put capacity study of MANETs, i.e. to derive the exact expressions for the 
throughput capacity of such networks. For the regular Manhattan and ring 
networks, Mergen and Tong [18] derived their throughput capacity in closed- 
form. Neely et al. [19] explored the exact capacity of cell-partitioned MA¬ 
NETs, and revealed a fundamental tradeoff between the throughput capacity 
and packet delay in such networks. Following this line, Wang et al. [20] ex¬ 
tended the tradeoff results in [19] to a multicast scenario. Gao et al. [21] 
considered a more general network scenario with a general setting of trans¬ 
mission range and derived its capacity under the group-based scheduling. 
Liu et al. further explored the exact throughput capacity of MANETs with 
packet redundancy [22] and power control [23]. Recently, Chen et al. studied 
the exact throughput capacity of MANETs with directional antennas [24] 
and explored the efficient approximations for the exact throughput capacity 
of MANETs with ALOHA protocol [25]. 

It is notable, however, that one common limitation of all these studies is 
that to make their analysis tractable, they all assume the relay buffer of a 
node, which is used for temporarily storing packets of other nodes, has an 
inhnite buffer size. This assumption does not hold for a practical MANET, 
where the buffer size of a mobile node is usually limited due to both its storage 
space limitation and computing capability limitation. Thus, for the practical 
capacity study of MANETs, the constraint on buffer size should be carefully 
addressed. Notice that the throughput capacity modeling with practical 
limited-buffer constraint still remains a technical challenge. This is mainly 
due to the lack of a general theoretical framework to efficiently characterize 
the highly complicated buffer occupancy behaviors in such networks. 

As a step to address above limitations, this paper studies the exact 
throughput capacity for a class of MANETs, where each node is associated 
with a shared and limited relay buffer to temporarily store the packets of 
other nodes [26, 27] and the flexible two-hop relay routing scheme is adopted 
for packet forwarding. The two-hop reply serves as an important routing 
protocol for practical MANETs [28], since it is simple and can be imple¬ 
mented easily in a distributed fashion. Also, the two-hop reply is efficient 
in the sense that it has the capability of achieving the throughput capac¬ 
ity under many important MANET scenarios [19, 21, 25]. For this class of 
buffer-limited MANETs, as a thorough extension of our previous work [29], 
this paper develops a general theoretical framework to enable the analytical 
study on their exact throughput capacity to be conducted. 
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The main contributions are summarized as follows: 


• For one buffer-limited MANET concerned in this paper, we first pro¬ 
vide theoretical analysis to reveal the inherent relationship between its 
throughput capacity and its relay buffer blocking probability (RBP). 

• For the analysis of RBP, a novel theoretical framework based on the 
Embedded Markov Chain Theory is then developed to capture the 
complicated queuing process in a relay buffer. With the help of the 
framework and also the Queueing Theory, the RBP under any exoge¬ 
nous input rate and the closed-form expression for exact throughput 
capacity are then derived. 

• Case studies are further provided under two typical transmission schedul¬ 
ing schemes to illustrate how our theoretical framework can be applied 
for exact throughput capacity analysis under one given transmission 
scheduling scheme. The corresponding capacity optimization issue and 
scaling law performance are also explored. 

• Finally, extensive simulation and numerical results are provided to 
demonstrate the efficiency of our theoretical framework on capturing 
the throughput behavior of a buffer-limited MANET and to illustrate 
the impacts brought by the buffer constraint. 

The remainder of this paper is outlined as follows. Section 2 introduces 
the system models, routing scheme and dehnitions involved in this study. We 
present the overall framework for throughput capacity analysis in Section 3, 
and then develop an Embedded Markov Chain-based framework in Section 4 
for the evaluation of RBP and exact throughput capacity. Section 5 deals 
with the case studies, capacity optimization issue and scaling law results, and 
Section 6 provides the simulation and numerical results. Finally, we provide 
the related work in Section 7 and conclude this paper in Section 8. 

2. Preliminaries 

This section presents the system models, buffer constraint, routing scheme 
and basic dehnitions involved in this study. Table 1 summarizes the main 
notations. 
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Table 1: Main notations 


Symbol 

Quantity 

n 

Number of nodes 

m 

The network is partitioned into m x m cells 

B 

Relay buffer size 

A 

Exogenous packet arrival rate 

PbiX) 

Relay-buffer blocking probability 

Te 

Throughput capacity 

a 

The probability that a node selects to execute S-R 
transmission when it gets access to the wireless channel and 
cannot execute the S-D transmission_ 


2.1. System Models 

As illustrated in Fig. 1 that we consider a time-slotted and cell-partitioned 
network with n mobile nodes and m x m non-overlapping equal cells, where 
the nodes roam from cell to cell over the network according to the indepen¬ 
dent and identically distributed (i.i.d) mobility model [11, 19], and during a 
time slot the total amount of data that can be transmitted from a node to 
another is fixed and normalized to one packet. With the i.i.d mobility model, 
each node independently selects a cell among all cells with equal probability 
at the beginning of each time slot and then stays in it during this time slot, 
so the location of each node is i.i.d and uniformly distributed over all cells in 
each time slot, and between time slots the distributions of nodes’ locations 
are independent. 

We adopt the i.i.d. mobility model here mainly due to the following rea¬ 
sons. First, the mathematical tractability of this model allows us to gain im¬ 
portant insights into the structure of throughput capacity analysis. Second, 
as illustrated in [19] and to be demonstrated in Section 6.2, the throughput 
capacity result derived under this model can be applied to other typical mo¬ 
bility models like the random walk model. Third, the analysis under this 
model provides a meaningful theoretical performance result in the limit of 
inhnite mobility [19]. Notice that the mobility model determines the distribu¬ 
tions of nodes’ locations and thus the opportunity that nodes encounter with 
each other, so it will affect the throughput capacity that can be achieved, as 
to be shown in Section 4 and Section 5. 

We consider the permutation traffic model widely used in previous liter- 
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Figure 1: System models. 

ature [12], where n distinct unicast traffic flows exist in the network, each 
node is the source of one flow and meanwhile the destination of another flow. 
We assume that the exogenous (self-generated) packet arrival at each node 
follows the i.i.d Bernoulli process with mean rate A packets/slot. 

2.2. Node Model and Buffer Constraint 

As illustrated in Fig. 1 that we consider a practical node model similar 
to that of previous studies on buffer-limited wireless networks [26, 27], where 
each node maintains two independent queues, one local queue with unlimited 
buffer size for storing the exogenous packets of its own flow and one shared 
relay queue with fixed size B for storing the relay packets coming from all 
other n — 2 traffic flows. The local queue follows the FIFO (flrst-in-flrst-out) 
discipline. The relay queue follows the quasi-FIFO (Q-FIFO) discipline, i.e., 
the packets destined to the same node in the relay queue follow the FIFO 
discipline. 

We adopt this node model here mainly due to the following reasons. First, 
in a practical network, each node usually reserves a much larger buffer space 
for storing its exogenous packets rather than the relay packets. Second, even 
though the local buffer space is not enough when bursty traffic comes, the 
upper layer (like transport layer) can execute congestion control to avoid the 
loss of local packets. Third, as illustrated in Fig. 1 that the local queue serves 
as an ingress for exogenous packets entering a MANET, so it is actually 
an external module of the MANET and its buffer size does not affect the 
maximal achievable throughput the MANET can support. Finally, the relay 
queue serves as an important internal module of a MANET for storing and 
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Figure 2: Illustration of the 2HR-a routing scheme. 

forwarding relay packets, so its buffer size has a critical impact on the network 
throughput performance. 

2.3. Routing Scheme 

To support the efficient operation of buffer-limited MANETs, we extend 
the Two Hop Relay (2HR) algorithm [6] and consider a more flexible Two Hop 
Relay scheme with parameter a (2HR-Q; for short), which incorporates both 
a control parameter a for transmission control and a handshaking mechanism 
for avoiding unnecessary packet loss. 

Regarding the 2HR scheme, for a tagged flow (S, D) with source node S 
and destination node D, when S gets access to the wireless channel in a time 
slot, it will transmit a packet directly to D (S-D transmission) if D is within 
its transmission range; otherwise with probability 0.5, S selects to transmit a 
self-generated packet to a relay node (S-R transmission), or deliver a packet 
of other nodes to the corresponding destination (R-D transmission). 

Notice that under the buffer-limited scenario, when node S executes the 
S-R transmission while the relay queue of the receiver is full, then the trans¬ 
mission will not be successful and the transmitted packet will be lost. To 
facilitate the operation of buffer-limited MANETs and improve the through¬ 
put performance, we adopt here the 2HR-Q; scheme, which is an extension of 
the 2HR scheme in the following two aspects. Eirst, as illustrated in Eig. 2 
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that we introduce a parameter a to flexibly control the probability that S 
selects to conduct S-R transmission, i.e., when S gets access to the wireless 
channel and node D is not within its transmission range, S selects to transmit 
a self-generated packet to a relay with probability a, and deliver a packet of 
other nodes to the corresponding destination with probability 1 — a. Thus, a 
represents the level of selfishness of a node, from 0 (fully selfless) to 1 (fully 
selfish), and it is expected that a should be set appropriately according to 
the network settings to achieve the optimal throughput performance. Sec¬ 
ond, to avoid the unnecessary packet loss in S-R transmission, the 2HR-Q; 
scheme further adopts a handshake mechanism to confirm the relay-buffer 
occupancy state of a receiver, where the S-R transmission will be conducted 
only when the relay queue of the intended receiver is not full. 

2.4- Basic Definitions 

Relay-buffer Blocking Probability (RBP); For a buffer-limited MANET 
with the 2HR-Q; scheme and a given exogenous packet arrival rate A to each 
node, the relay-buffer blocking probability of a node is defined as the 

probability that the relay queue of this node is full. 

Throughput: The throughput of a flow is defined as the time average of 
number of packets that can be delivered from its source to its destination. 

Throughput Capacity and Optimal Throughput Capacity; For a 
buffer-limited MANET with the 2HR-Q; scheme, its throughput capacity 
is defined as the maximal achievable per-flow throughput the network can 
support, independent of the arrival rate. The optimal throughput capacity 
is defined as the maximum value of throughput capacity Tc optimized over 
the control parameter a, i.e., T* = max T^. 

aG[0,l] 

3. Throughput Capacity Analysis 

Eor a buffer-limited MANET with the 2HR-Q; scheme, we denote by psd, 

Psr and Prd the probabilities that in a time slot a node gets access to the 
wireless channel and selects to execute S-D, S-R and R-D transmission re¬ 
spectively. With the help of these basic probabilities, we can establish the 
following theorem regarding the throughput capacity of the network. 

Theorem 1. For a buffer-limited MANET with the 2HR-a scheme, its through¬ 
put capacity Tc is determined as 


Tc =Psd + Psr{^ -PbW) packets/slot, 


( 1 ) 



where X is the unique solution of the following equation 

X=Psd+Psr{^-Pb{X)). (2) 

Proof. To prove the Theorem, we hrst demonstrate that there exists an 
unique solution A for the equation (2), and then show that the through¬ 
put is A when packet arrival rate A < A, but the throughput is always A 
when A > A. 

From our system models, it is clear that each flow experiences the same 
service process without priority, so the queuing process of each flow is iden¬ 
tical and we can focus on a tagged flow in our analysis. Under the 2HR-Q; 
routing scheme, the delivering process of a packet from its source to destina¬ 
tion involves at most two stages. The first stage is the queuing process at its 
source node, while the second stage is the queuing process at one relay node 
if the packet is not directly delivered to its destination. 

Concerning the first stage queuing process, the local queue there can be 
modeled as a Bernoulli/Bernoulli queue with arrival rate A and service rate 
p.s{X) determined as 

Ps{X) = Psd +PsrO-- Pb{X)) . (3) 


We can easily see that: 1) when A = 0, we have P6(0) = 0, so Ps(0) = 
Psd + Psr > A; 2) as A increases, Pb{X) tends to increase, leading to a decrease 
in Ps(A); 3) when A = Psd + Psd we have Pb{X) > 0, p.s(A) < A. Based on 
these properties of service rate Ps{X), we know that there exists an unique 
0 < A < Psd + Psr such that A = Ps{X). 

Considering a time interval [0, t], we denote by mo(t) and mi(t) the num¬ 
ber of packets being buffered in all local queues and all relay queues at time 
slot t, respectively. Since the total number of exogenous arrival packets dur¬ 
ing this interval is nXt, then the throughput Th is determined as 


^ nXt - mojt) - mijt) 

t^Qo n ■ t 


( 4 ) 


Since the relay-buffer of each node has a hxed size B, then ttHfH < q 

lim = 0. 

t^oo 

For the case A < A, we denote by Lg the queue length of the local queue, 
then the expectation EjL^} of Lg is given by [30] 


E{Lg} 


A-A^ 

psi,X^ 


( 5 ) 
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Since when A < A, we have yUs(A) > A, so the queue length is bounded 

in this case. Thus, lim = 0 and Th = A. 

When A > A, then /is(A) < A, leading to an increasing number of packets 
buffered in the local queues. By applying the law of large numbers, we have 
that as f —)■ cxD 

^o{t) a.s 

— -^ n(A -/i,(A)). 

Based on (4), we then have Th = X when A > A. 

Thus, the throughput capacity of the concerned network is determined 
as 

Tq hsi^X^ Psd T Psr Pb(^X)''j . 


□ 


Remark 1. Notice that for the heterogeneous network scenario, the network 
level capacity region is dehned by a vector A = {Ai, A 2 , • • • , A„}, where Aj 
denotes the maximum packet arrival rate of node i that the network can 
stably support. For the homogeneous network scenario considered in this 
paper, the behavior of each node is the same, so we have Ai = A 2 = ■ ■ ■ = A„. 
Thus, the network level capacity here reduces to the per node capacity as 
expressed in (1). 


4. Embedded Markov Chain Framework 

The result in Theorem 1 indicates that for the throughput capacity anal¬ 
ysis of the concerned MANET, we need to determine RBP pb{X) in such net¬ 
work. To address this issue, in this section we hrst utilize a two-dimensional 
Markov Chain (MC) to depict the complicated state transitions of a relay 
queue, then convert the two-dimensional MC into a new Embedded Markov 
Chain (EMC) to obtain its one-step transition probability, such that a com¬ 
plete theoretical framework is developed to enable the RBP and thus the 
exact throughput capacity of the concerned MANET to be derived. 

4.T State Machine of Relay Queue 

Without loss of generality, we focus on the relay queue of a general node 
and use a two-tuple (i, k) to denote the state that the relay queue contains i 
packets and these packets are destined for k distinct destination nodes, here 
0<i<i?, A; = 0 when i = 0 and 1 < k < i when i > 0. As illustrated in 
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(b) State machine of a relay queue. 

Figure 3: Transitions and state machine of a relay queue. 


Fig. 3a that when the relay queue is in state (i, k) at current time slot, only 
one of the following transitions may happen in the next time slot: 

• Case 1: transition from {i,k) to {i + 1,/t), i.e., a packet enters the 
relay queue and this packet is destined for a destination same as one 
of packet (s) already in the queue. 

• Case 2: transition from {i,k) to {i + l,k + 1), i.e., a packet enters 
the relay queue and the destination of this packet is different from all 
packet (s) in the queue. 

• Case 3: transition from (i, k) to {i — 1, A;), i.e., a packet from the relay 
queue is delivered to its destination, but there still exist other packet (s) 
in the relay queue destined for this destination. 
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Figure 4: State machine of the EMC. 

• Case 4: transition from [i, k) to {i — l,k — 1), i.e., a packet from the 
relay queue is delivered to its destination, and none of the remaining 
packet (s) in the relay queue is destined for this destination. 

• Case 5: transition from {i,k) to {i,k), i.e., no packet entering into or 
departing from the relay queue. 

Based on above transitions, the state machine of the relay queue can then 
be modeled as a discrete-time hnite-state Markov chain illustrated in Fig. 3b, 
where flj denotes the set of states {i,k),k < i. 

4-2. Markov Chain Collapsing 

It is easy to see that the total number of states of the Markov chain in 
Fig. 3b is 1 -|- , so the Markov chain will become too complicated to 

be solved when B is big. Notice also that for the analysis of RBP ^^(A), 
the limiting distribution on each set j rather than the limiting distribution 
on each state {i, k) is of concern. For these reasons, we apply the novel 
“Markov Chain Collapsing” [31] technique to convert the Markov chain in 
Fig. 3b into a new EMC illustrated in Fig. 4, where flj denotes the state that 
the relay queue contains i packets and denotes the one-step transition 

probability from flj to Qj. 

With the help of the EMC model in Fig. 4 and also the Markov Chain 
model in Fig. 3b, we can establish the following lemmas regarding the eval¬ 
uation of and also the limiting distribution of the EMC. 

Lemma 1. The one-step transition probability of the EMC in Fig. 4 
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is determined as 


f Ps{^) ■ Psri 

i 




= < 


Q I • ■ Prd-} 
n — 6 + t 

1 Pfli,fli+i PQi,fli-i} 

0 , 


j = i + I < B 


j = i — 1 > 0 

j = i 
others 


where p.(\) = = JM' 

Proof. See Appendix A for the proof. 


( 6 ) 


□ 


Lemma 2. The limiting distribution 11 = {7r(r2o), of the 

EMC exists and is unique, and it is given by 


7r(fio) 




1 

Ef=oC.-/9*-Ps(A)*’ 

Q-/ 3 ^-p,{Xy 

j:f=oCh-/3^-Ps{Xy' 


0<i< B 


(7) 

( 8 ) 


where Q = (” and /3 = ^ = 

^ \ t / ^ Prd I—a 

Proof. See Appendix B for the proof. □ 


4 . 3 . Derivation of pb{X) and Tc 

Based on above EMC-based framework, we now provide analysis on the 
RBP Pb{X) and the exact throughput capacity T^, as summarized in following 
theorem. 


Theorem 2. (Main result) For a coneerned MANET with n mobile nodes, 
where each node is allocated with a relay-buffer of fixed size B and the 2HR-a 
scheme is adopted for packet delivery, its RBP is determined by the following 
equation 

Gb ■ ■ Ps{X)B 

Ef=oC*-/5*-p.(A)*’ 

and its throughput capacity Tc is determined as 

= + ( 10 ) 
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Proof. Notice that when the EMC is on state Q.b, the relay queue is full. 
It indicates that RBP Pb{\) is just equal to the limiting distribution 7r(r2B), 
so (9) follows. Since Ps(A) = ^ ’ equation (9) contains only one 

unknown quantity Pfe(A). By solving equation (9), we can then obtain the 
RBP Pb{\) under any exogenous packet arrival rate A. 

From Theorem 1 we know that as A approaches A, Ps(A) tends to 1. 
Substituting Ps(A) = 1 into (9), we have 


Pfe(A) 


Cb-P^ 

Ef=oC.-/3*' 


The formula (10) then follows by substituting (11) into (1). 


( 11 ) 


□ 


Remark 2. It is notable that our EMC-based framework for throughput 
capacity analysis is very general in the sense that it can be applied to other 
network models (like the continuous network model [25]) and other mobility 
models (like the random walk mobility model), as long as they lead to the 
same steady distributions of the nodes’ locations as the i.i.d mobility model. 


Based on Theorem 2, we have the following corollaries (see Appendix C 
for the proofs). 

Corollary 1. For a network with n >3, its throughput capacity Tc increases 
as relay-buffer size B grows. 

Corollary 2. With the setting of a = 0.5, i.e., each node executes S-R and 
R-D transmission with egual probability, Tc is determined as 

B 

Tc = Psd+Psr -(12) 

n — 2 + B 

Corollary 3. When the relay-buffer size B tends to infinity, the throughput 
capacity Tc is determined as 


'^c I B^oo 


Psd T Psr: 
Psd T Prd: 


a < 0.5 
a > 0.5 


(13) 
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5. Case Studies 

The results in Theorem 2 indicate for the evaluation of throughput ca¬ 
pacity of a MANET concerned in this paper, we need to determine the prob¬ 
abilities Psd and Psr in the network, which are further related to the schemes 
adopted for transmission scheduling. To demonstrate the applicability of our 
Embedded Markov Chain-based framework for throughput capacity analysis, 
this section provides case studies under two typical transmission scheduling 
schemes, i.e., local transmission scheduling (LTS) and group-based trans¬ 
mission scheduling (GTS). The related issues of capacity optimization and 
scaling law performance will be also explored. 

5.1. Throughput Capacity under LTS 

We hrst provide analysis of throughput capacity under the LTS scheme 
[19], where two nodes within the same cell can forward one packet during 
a time slot, and nodes within different cells cannot communicate. Same as 
[19], the node/cell density d is assumed to be 0(1) (independent of n) and 
without loss of generality, n is assumed to be even and source-destination 
pairs are composed as: 1 -H- 2, 3 -H- 4, • • • , n — 1 -H- n. 

We denote by po and pi the probabilities that there are at least two nodes 
in a cell and there is at least one source-destination pair in a cell, respectively. 
Based on the results of [19], we then have 



At a time slot, the total transmission opportunity in the network is -po, 
which is shared equally by all nodes, so we have 


n ■ {psd + Psr + Prd) = ■ Po- 


Similarly, 


n-Psd ^ -Pi- 


(16) 

( 17 ) 
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Combining with — = we have 

Prd J- ^ 


Psd = ^Pl, 


a 


Psr = -^{PO-Pl), 


a 


Prd 


d 


(po-Pi)- 


( 18 ) 

(19) 

( 20 ) 


Substituting (18) and (19) into (10), we can see that the throughput 
capacity under the LTS scheme is determined as 


a 


Tc = -:Pi + -jipo-pi) 1 - 




B 


d 


d 




( 21 ) 


Remark 3. From Corollary 3 we can see that when a = 0.5 and B —)■ oo, 
then (21) is reduced to the capacity result in [19], i.e., Tc = 


Regarding the optimal throughput capacity T* and the corresponding 
optimal setting of a* under the LTS scheme, we have the following theorem. 

Theorem 3. For a concerned MANET with the 2HR-a scheme and a fixed 
relay-buffer size B, its optimal throughput capacity T* under the LTS scheme 
is determined as 




1 po - Pi Ki*) 
d^^ d{l + M)h{T) T Cb' 


( 22 ) 


and the corresponding optimal transmission ratio a* is given by a* = 
where 

HX = E C. ■ 9®". (23) 

i=0 

h'ipf) is the derivative ofhp), andr* is determined by the following equation 


h{l*)[h{P) + Cs] = (1 + P)CBh\P). (24) 

Proof. We dehne 7 = 1^ (i.e., « = ^, /3 = i), and ^( 7 ) = (1 + 
7 ) ^1 + From (21) we can see that the optimal throughput capacity 
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T* is determined as 


rjn* 


max Tc 
aG[0,l] 


-Pi + -:{po-Pi) ■ max 
a a o 6 [o,i] 

1 Po-Pi 1 
« d min ^( 7 )’ 

7>0 



1 - 


Cb-P‘ 


Ef.o c. ■ P‘ 


(25) 


We can see that: 1) g{'y) is an Elementary Function [32], so it is continuous 
and differentiable on the interval 7 G ( 0 , cxd); 2 ) lim 5 f( 7 ) —)■ cxd, lim g{^) 

7^0 ■y—yoo 

00 and gl'y) > 0. According to the Extreme Value Theorem [33], there must 
exists 0 < 7 * < cxD such that 0 < ( 7 ( 7 *) < g{'y) for V 7 G (0, cxd) and g'i'y*) = 0, 
so equation (24) follows. Then formula (22) follows by substituting 7 * into 
(25). 

□ 


Based on Theorem 3 we have the following corollary. 

Corollary 4. For any setting of n and B, a* < 0.5; when B —)■ oc, 
= 0.5 and 

Proof. See Appendix D for the proof. □ 

With the help of exact expression of throughput capacity (21), the achiev¬ 
able scaling law under the LTS scheme can be further explored, as shown in 
the following corollary. 

Corollary 5. When the number of nodes in the concerned MANET tends to 
infinity, each node can achieve a throughput Tfin) as 

TPn) = e ([®]) . (26) 

Proof. See Appendix E for the proof. □ 

Remark 4. Notice that Corollary 5 is the first time to reveal an achievable 
order sense (in 0 form) on per node throughput of a buffer-limited MANET, 
not just an upper bound (in O form) provided in [26], and it indicates that as 
n increases, the concerned MANET can still achieve a non-vanishing through¬ 
put as long as its relay-buffer size B grows at least linearly with n. 
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5.2. Throughput Capacity under GTS 

We further conduct analysis on throughput capacity under the GTS 
scheme [12, 21, 34], With GTS scheme, all cells are divided into different 
groups, where any two cells in the same group have a horizontal and vertical 
distance of some multiple of e cells. Thus, the MANET has groups and 
each group contains J = \wS/e^\ cells. Each group becomes active every 
time slots and a node in an active cell can transmit one packet to an¬ 
other node within a horizontal and vertical distance oi v — 1 cells. Based 
on the Protocol Model [5], to ensnre that concnrrent transmissions are not 
interfering with each other, e is determined as [ 21 ] 

e = min{|’(l-I-A)\/2z/-I-(27) 


Gonsidering a given time slot and a given active cell c, we denote by ps 
the probability that there are at least one node within c and another node 
within the transmission range of c, and denote by p^ the probability that 
there are at least one sonrce-destination pair within the transmission range 
of c and for each of such pair(s), at least one of its two nodes is within c. 
Based on the results of [21], we have 

-1)" - (28) 
TDj 

p- = Arl™"" - - 2' + 1)"'% (29) 

m 

where / = (2z/ — 1)^. Notice also that psd = ^^ 4 , p^r = — Pa) and 

Prd = — P 4 ). By snbstitnting these results into (10), the thronghpnt 

capacity under the GTS scheme is then determined as 


Tc 


■J aJ 

-Pa H-(P3 

n n 


Pa) 



Ef=„c.-/3‘) 


(30) 


We can see that when a = 0.5 and B ^ 00 , then (30) is rednced to the 
capacity resnlt in [21], i.e., Tc = Based on the proofs similar to that 

of Theorem 3 and Gorollary 5, we have the following corollary regarding the 
optimal thronghpnt capacity T* and scaling law nnder the GTS scheme. 


Corollary 6. For a concerned MANET, its optimal throughput capacity T* 
under the GTS scheme is determined as 


JjPS-PA) h(7*) 
n{l + M) hlM) + Cb^ 


(31) 
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and the corresponding optimal transmission ratio a* is given by a* = 
where h{'y) and 7 * are determined by (23) and (24), respectively. The order 
of throughput capacity Tdn) under the GTS scheme is determined as 

T,(n) = 0 (f) • (32) 

6. Simulation Results and Discussions 

In this section, we first provide the simulation results to validate our 
theoretical framework for the throughput capacity analysis of buffer-limited 
MANETs, and then apply our theoretical results to illustrate the performance 
of such networks. 

6.1. Simulation Setting 

For the validation of our framework, a C-I--I- simulator was developed 
to simulate the packet delivery process in the concerned MANETs under 
both LTS and GTS schemes [35]. In addition to the i.i.d mobility model, 
the random walk model was also implemented in the simulator. Under the 
random walk model, at the beginning of each time slot, every node indepen¬ 
dently selects a cell among its current cell and its 8 adjacent cells with equal 
probability 1/9, then stays in it until the end of this time slot [ 8 ]. 

Two network scenarios of (n = 72, m = 6, B = 5,a = 0.5) and (n = 
200, m = 10,5 = 8 , 0 ; = 0.3) are considered in the simulation, where we set 
u = 1 and A = 1 for the GTS scheme [36]. To simulate the throughput, 
we focus on a specihc node and count its received packets over a period of 
2 X 10® time slots, and then calculate the averaged number of packets this 
node can receive per time slot. The system load p is dehned as p = A/Tc, 
and Tc is given by (21) and (30) for the LTS and GTS, respectively. 

6.2. Simulation Results 

To validate the throughput capacity results (21) and (30), we provide 
plots of throughput versus system load p in Fig. 5. It can be observed from 
Fig. 5 that the simulation results agree well with the theoretical ones under 
both LTS and GTS schemes, indicating that our framework is highly efficient 
in capturing the throughput behaviors of concerned buffer-limited MANETs. 
We can see from Fig. 5 that just as Theorem 1 predicates that for a concerned 
MANETs, its throughput increases linearly with p when p < 1 and then keeps 


19 




(a) Throughput performance under the LTS 
scheme. 



(b) Throughput performance under the GTS 
scheme. 


Figure 5: Throughput performance under the LTS and GTS schemes. Case 1: n = 72, m = 
6, B = 5, a = 0.5. Case 2: n = 200, m = 10, i? = 8, a = 0.3. 

as a constant T^. determined by (1) when p > 1. A further observation of 
Fig. 5 indicates that for a network under the random walk mobility model, 
its throughput performance is very similar to that under the i.i.d mobility 
model. This is expected since according to Remark 2 and [19, 37], nodes in 
the network with i.i.d mobility model and random walk model have the same 
steady distribution, leading to the same throughput capacity performance 
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0.16 



Figure 6: Throughput capacity Tc vs. relay-buffer size B. 



Figure 7: Throughput capacity Tc vs. number of nodes n. 

under these two mobility models. 

6.3. Throughput Capacity 

With the help of our theoretical results, we illustrate here the impacts of 
network parameters on the throughput capacity. Notice that for a concerned 
MANET its overall throughput behavior under the LTS is very similar to that 
under the GTS, so we consider the LTS scheme here only for illustration. 

We hrst summarize in Fig. 6 how throughput capacity Tc varies with 
relay-buffer size B under two network scenarios of (n = 72, m = 6 ) and 
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Figure 8: Throughput capacity Tc vs. transmission ratio a. 

{n = 200,m = 10), where a is fixed as 0.5. We can see from Fig. 6 that 
just as discussed in Corollary 1, the throughput capacity of a buffer-limited 
MANET can be improved by adopting a larger relay-buffer in such network. 
A careful observation of Fig. 6 shows that as B increases the capacity Tc first 
increases quickly and then gradually converges to a constant determined by 
Corollary 3. This observation indicates we can determine a suitable buffer 
size B according to the requirement on network capacity such that a graceful 
trade-off between capacity performance and buffer cost can be achieved. 

To further illustrate the impact of buffer constraint on the throughput 
capacity, we show in Fig. 7 the relationship between Tc and n under three 
typical relay-buffer settings, i.e., B is fixed as a constant (5 here), B = n 
and B ^ oo. We can see from Fig. 7 that in general Tc decreases as n 
increases, but as n —)■ cxd, Tc vanishes to 0 when B is fixed, while a non-zero 
constant throughput capacity can still be achieved when B = n or B ^ oo. 
These behaviors are expected, since the result in Corollary 5 indicates that 
to achieve a non-vanishing throughput capacity in the concerned MANET, 
its relay-buffer size B should grow at least linearly with n. 

6.4- Optimal Throughput Capacity 

To illustrate the optimal throughput capacity, we show in Fig. 8 the 
impact of transmission ratio a on throughput capacity Tc under the settings 
of n = 72, m = 6 and B = {5,20}. We can see from Fig. 8 that under a 
given setting of T, as a increases Tc first increases and then decreases, and 
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(a) Optimal transmission ratio a* vs. relay-buffer 
size B. 



(b) Optimal transmission ratio a* vs. number of 
nodes n. 


Figure 9: Optimal transmission ratio a* vs. B and n. 

just as discussed in Theorem 3 that there exists an optimal a* to achieve 
the optimal throughput capacity T* . This is mainly due to the reason that 
the effects of a on are two folds. On one hand, a larger a will lead to 
a higher probability of conducting S-R transmission; on the other hand, a 
larger a will result in a higher RBP thus a lower opportunity of conducting 
the S-R transmission. As a snmmary, in order to improve the throughput 
performance of a bnffer-limited MANET, nodes shonld cooperate with each 
other, and they shonld be neither too selfish nor too selfless. 
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Based on the results of Theorem 3, we illustrate in Fig. 9 how the optimal 
transmission ratio a* is related to B and n. We can see that just as proved in 
Corollary 4 that a* increases as B grows while it decreases as n grows, and the 
optimal transmission ratio never exceeds 0.5. These behaviors indicate that 
in a network with the hxed number of nodes n, if we upgrade the capacity of 
each node by adopting a larger relay-buffer, we should accordingly allocate a 
higher probability for S-R transmission (i.e., nodes should be more selhsh), to 
achieve the optimal throughput capacity. On the other hand, when the relay- 
buffer size of each node is fixed, if we increase the scale of the network by 
accommodating more nodes, we should accordingly increase the probability 
of R-D transmission (i.e., nodes should be more selfless), to release the relay- 
buffer space and thus guarantee the optimal throughput capacity there. 

7. Related Work 

A signihcant amount of work has been devoted to the performance anal¬ 
ysis of wireless ad hoc networks, among which some initial studies focused 
on the buffer-limited network scenarios. 

Herdtner and Chong [26] explored the throughput-storage tradeoff in MA- 
NETs and showed that the per node capacity under the hnite buffer con¬ 
straint cannot achieve 0(1) even though node mobility is utilized. Later, 
Eun and Wang [38] considered TCP/AQM system and developed a doubly- 
stochastic analytical model to study the tradeoff among link utilization, 
packet loss, and buffer size. Subramanian et al. [39] developed a theoreti¬ 
cal framework for throughput analysis in Delay Tolerant Networks (DTNs). 
They derived closed-form expressions for the per node throughput capacity 
of a sparse network which consists of one source-destination pair and several 
mobile relay nodes with limited buffer size, and further extended their re¬ 
sults to the multi-cast network scenarios [40]. Wang et al. [41] studied the 
node buffer occupancy behaviors in static random wireless networks with 
intermittent connectivity and provided some scaling results of fundamental 
achievable lower bound for the occupied buffer size. Krifa et al. [42] fo¬ 
cused on the buffer management policies in DTNs. They demonstrated that 
drop-tail and drop-front policies are sub-optimal and meanwhile proposed 
an optimal buffer management policy based on global knowledge about the 
network. Following this line, Elwhishi et al. [43] developed a new message 
scheduling framework for Epidemic [44] and Spray-Wait [45] forwarding rout¬ 
ings in DTNs with hnite buffer to optimize either the message delivery ratio 


24 



or message delivery delay. 


8. Conclusion 

In this paper, we first revealed the inherent relationship between the 
throughput capacity and relay-buffer blocking probability in a buffer-limited 
MANET with the 2HR-Q; scheme, and then developed an Embedded Markov 
Chain-based framework to fully characterize the complicated packet delivery 
processes of the concerned MANET. Based on this framework, we derived the 
throughput capacity in closed form and further conduct cases studies under 
two typical transmission scheduling schemes to illustrate the impacts of some 
key network parameters on throughput capacity. It is expected the theoret¬ 
ical framework developed in this paper will be also helpful for exploring the 
throughput capacity of buffer-limited MANETs under other mobility mod¬ 
els and other transmission schemes. The results in this paper indicate that 
in large scale MANETs a non-zero constant throughput capacity can still be 
guaranteed as long as its relay-buffer size grows at least linearly with network 
size. Another interesting Ending of this paper is that for throughput capacity 
optimization in such MANETs, the optimal setting of transmission ratio in 
the 2IIR-Q; scheme there increases with the relay-bufier size but decreases 
with the network size, and it never exceeds 0.5. 

Notice that the theoretical framework and closed-form results for per node 
throughput capacity developed in this paper is based on the i.i.d mobility 
model, so one of our future research directions is to develop theoretical models 
for other more realistic mobility models, like the inter-meeting time based 
mobility model. 

Appendix A. Proof of Lemma 1 

Notice that each local queue is a Bernoulli/Bernoulli queue with exoge¬ 
nous packet arrival rate A, so the output process of the local queue is also a 
Bernoulli flow with rate A due to the reversible property of Bernoulli/Bernoulli 
queue [30]. 

Based on the property of the i.i.d mobility model we know that for a 
specific node, except itself and its destination, each of the remaining n — 2 
nodes will equal likely to serve as its relay. Similarly, for each node serving 
as a relay, except itself and its source, all the remaining n — 2 nodes will 
equal likely to forward packets to it. From (3) we know that the ratio of 
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S-R transmission to the service rate /is(A) is Hence, the packet 

arrival rate at a relay qnene Xr is determined as 


Ar = n - 2 A-—-- [n - 2 

hs(A) 

= Ps(A)Psr (1-Pb(A)). (A.l) 

We denote by P(i,k),ni+i the transition probability from state {i, k) to set 
Hj+i, 0 < i < -B. Since when a relay qnene is full, its packet arrival rate will 
be 0. Thus, we have 


V{i,k),Qi+i ■ (1 - Pb(A)) + 0 ■ Pfe(A) = Xr, 

=^P(i,fc),Oi+i = Ps(A) ' Psr- (A-2) 

Notice that the transition probability pni,Qj of the EMC is the “set- 
averaged” transition probability of the original Markov chain [31], then we 
have 


Pni,nj = ^P{i,k),Qj ■ P ((f,/c)lfli), (A.3) 

k=l 

where P ((f, k)\kli) is the conditional probability that relay queue is in state 
(f, k) given that it belongs to the set fl*. Substituting (A.2) into (A.3) we 
have 

PUi,Qi+l PsiX) ■ Psr- 

A node gets a R-D transmission opportunity with probability prd, and 
due to the i.i.d mobility model, this opportunity arises for each of the n — 2 
destination nodes with equal probability. Thus, we have 

= k • (A.4) 

Based on (A.3), it is clear that in order to obtain pni,ni.i, we should derive 
the conditional probability P ((f, /c)|r2j). To address this issue, we utilize the 
Occupancy approach (Chapter 1 in [46]). When the relay queue has i packets 
being buffered, where each packet may be destined for any one of the n — 2 
destination nodes, the number of all possible cases is given by 
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- 


(A.5) 


- 3 + r 
i 

Suppose that these i packets are destined for k distinct destination nodes, 
the number of cases is given by 


E{i^k) 


n — 2 
k 


(k-l) + {i- k) 
i — k 


(A.6) 


Due to the i.i.d mobility model, each of these cases occurs with equal proba¬ 
bility. According to the Classical Probability, P ((i, /i:)|Dj) is then determined 
as 





(A.7) 


It can be easily verihed that P [{i,k)\kli) = 1. Substituting (A.7) into 

k<i 

(A.3) we have 


— 1 


i 




• (fc-1) kprd 1 
(n-3+i^ ‘ n - 2 J 





i 

n — 3 + i 


■ Prd- 


Appendix B. Proof of Lemma 2 

Based on the Lemma 1, the one-step transition probability matrix P of 
the EMC is given by 


Poi,no Poi,Oi Veil,0.2 


(B.l) 


PQ.b->^B — 1 P^b->^b 
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From the transition matrix (B.l) and the state machine in Fig. 4, we can 
see that: 1) the EMC is irreducible', 2) each state 12* is recurrent, 0 < i < B', 
3) the period of each state is 1, so each state is aperiodic. Based on these 
properties, we know that the EMC is ergodic, so its limiting distribntion 
exists and is nnique, which is jnst the stationary distribntion of the EMC 
[47]. Thus, we have 


n p = n, 

PsW ■ Psr ■ 7r(fio) = ■ vr(fii), 

n — 2 

Ps(A) ■ Psr ■ 7r{Qi) = ^ . 7r(n2), 

^ n-1 

* * * ? 

Ps{X) ■ Psr ■ ^ • 7r(fiB), 

n — 6 + B 

=s-!r(nj) = Cj ■ ff ■ ■ !r(no), 


(B.2) 


(B.3) 


Combining (B.3) with the normalization equation 
results (7) and (8) then follow. 


Appendix C. Proofs of Corollaries 1, 2 and 3 
Proof of Corollary 1: Let ai , then 







ELo - 2 +/c)Ci • 


1 + k + EiLo + l)Ci+i • 

Since 

71 — 2 -1- ? 

(k + l)Ci+i — (/c + 1) • Q, 

% +1 

and 

(/c + l)(n — 2 + i) — (i + l)(n — 2 + /c) 
= {n — 3){k — i) > 0, 
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then 


{k + l)Cj+i > (n — 2 + k)Ci, 


Substituting the result into (10), the Corollary 1 then follows. 

Proof of Corollary 2: When a = 0.5, then 13 = 1, and (10) is simplihed 
as 


Since 


T 

C 


Psd Psr 



(C.l) 


S [(« - 3) X (n - 4) ■ ■ • X 1 

+ (n - 2) X • • • X 2 + ■ ■ ■ + (n - 3 + 5) X • • ■ X (5 + 1 )] 
1 {n — 2 + B) X ■ ■ ■ X {B + 1) 

(n — 3)! n — 2 

/n — 2 + B\ 

= [ B 


(C.2) 


substituting (C.2) into (C.l), the Corollary 2 then follows. 

Proof of Corollary 3: For the case a = 0.5, since lim —= 1, 

substituting it into (12) we have 


Tc = Psd + Ps 
a=0.5,B^oo 


For the case a < 0.5, we have /3 < 1 and 

1 


B 

E 

i=0 


Q • = 


(n-3)! 

X [(n — 3) X • • • X 1 X + (n — 2) X • ■ • X 2 X 
H- \- {n-3 + B) X ■ ■ ■ X (5 + 1) X f3^] 

■ + / 3 ^-^ + • • • + 


(n-3)! 
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(n-Z+B \ 

(§') 

_ j^n— 2 +B \ 3) 

^ (n-3)! V 1-/9 ) 

where /(/9)^^^ denotes the k-th order derivative of f{f3). Since 

lim 1 - = 1, 

B^oo 


we have 


B 1 / 1 \ 

lim V Ci • = --- -- 

“ 3)! Vl - 
1 


(1-/9) 

and then 


lim Cs/3^(1< lim {B + nTl3^ 

B^oo B^oo 

< lim 2^B^I3^. 

B—^oo 


Since 

79! 

lim = lim ^ = lim ' = 0, 

tc—>-00 tc—>-00 ^ tc—>-00 — In py^ ‘ ^ 

substituting it into (10) we have 


Tq Psd Psr 

a<0.b,B^oo 


For the case a > 0.5, we have (3 > 1 and 

^ Cb • EfJo c, • 

Z)^0 Ci • /5* 1 + /3 Cj+1 • /3* 

1 


1 _l_ a Ci+i-/?^ 

Ef^V Ci-/3» ^ ^ Ef=“o' Ci-(Z* 


(C.3) 


(C.4) 


(C.5) 
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Since 


then 


Ef=-o' c. • 

lim = 1 

Ef=o'C.-/3* 


, , Cb-/3^ 1 

/3>l,B^oo Ef=oC*-/9* ^ 


Substituting it into (10) we have 


(C. 6 ) 


Tc Psd “1“ Psr n 

a>0.5,-B—)-oo P 

= Psd + Ps 


a 


1 — q; 


— Psd “1“ Prd- 


Appendix D. Proof of Corollary 4 

Considering 7 G (0,1], the first order derivative of g{;^) is 


+ Cb] - (1 + p)CBh'{-i)}. 

h[pY ■' -' 

(a) 

For Vn > 3, when i? = 1, (a) is determined as 

(a) = 7(7 + n - 2) - (1 + 7 )(n - 2) = ( 7 ^ - 1) - (n - 3) < 0. 
When B = k, we assume that 

(a) = hk{hk + Cfc) — (1 + p)Ckh'j^ < 0 , 
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where hk and h'f. are the abbreviations of and h'{'y) under B = k, 
respectively. When B = k + 1, we have 


Since 


(a) = hk+i{hk+i + Cfc+i) — (1 + 'y)Ck+ih'j^^i 
= 7 ■ i^k + Cfc) • [y{hk + Cfc) + Cfc_|_i] 

— (1 + 7 ) ■ Cfc+i • [hk + 'jh'f, + Cfc] 

= 'j’^hkjhk + Ck)^ +'y'^Ck{hk + Ck) + yCk+ikk^ 
(h) (ci) 

+ jCfcCfc+i — (1 + 7)Cfc-|-ihfc — 7(1 4 - '^)Ck+ih'k 

{dl) (C 2 ) ( 62 ) 

— (1 + 7 )CfcCfc+i. 

'-V-' 

id2) 


(bi) - (62) = 7 [ 7 ^fc(hfc + Ck) - (1 + 7)Cfc+ihy 
< 7[/ifc(/ifc + Ck) — (1 + 'y)Ckh'j^] < 0, 

combining (ci),(c 2 ) and {di),{d 2 ) we have 

(a) < 7 ^Cfc(/ifc + Ck) — Ck+ihk — CkCk+i 

= [hk + Cfc)( 7 ^Cfc — Cfc+i) < 0 . 

According to the above mathematical induction, we can conclude that g'{'y) < 
0 for 7 G ( 0 , 1 ) and g'iX) < 0 . Thus, gipf) monotonically decreases when 
7 G (0,1], so we know that 7 * > 1 and a* = < 0.5. 

For the limiting case B —)■ 00 , from (13), (19) and (20) we can easily see 
that a*\B^oo = 0.5 and T^*|s^oo = 

Appendix E. Proof of Corollary 5 

As n ^ cxD we have 


Po = 1 - 

\ 

^ 1 - e”"* - de”"*. 


1-1 

n 


\-d 


1-1 

n 


\-d-l 
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Pi = 1 - 


1 - 


dP 




= 1 - e2' 




^ 1 - e° = 0, 


Cb ■ Efco* Cj ■ /?• 

Ef.„Q./3- Ef=„c.'/3' 


-)■ 


B 


(n — 3 + -B)/9 

Substituting above results into (21), we have that as n —)■ cxd 




q;(1 




de 


-d\ 


B 


djS 


(n — 3 + -B) 


= 0 


5 


n 
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